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We show how to perform integrals over products of distri- 
butions in coordinate space such as to reproduce the results 
of momentum space Feynman integrals in dimensional regu- 
larization. This ensures the invariance of path integrals under 
coordinate transformations. The integrals are performed by 
expressing the propagators in 1 - e dimensions in terms of 
modified Bessel functions. 



I. INTRODUCTION 

While the Schrodinger equation and time-sliced path 
integrals, as defined in the textbook |Q, are invariant 
under coordinate transformations by construction, this 
invariance is completely nontrivial in perturbatively de- 
fined path integrals in curvilinear coordinates. The rea- 
son for this is that the Feynman integrals appearing in 
the perturbation expansion are multiple temporal inte- 
grals over highly singular integrands, which are products 
of generalized functions [for general discussion, see || ]. 
In the presently existing theory of distributions, gener- 
alized functions can only be combined linearly and by 
a convolution, but products of them are still forbidden 
calling for regularization ||. Thus, if we want to de- 
scribe quantum mechanics in curvilinear coordinates by 
perturbative defined path integrals, we must extend the 
theory of distributions to allow for integrals over prod- 
ucts of generalized functions. This will extend the linear 
space of distributions to a group. 

In a previous note pj we have shown that invariance 
can be achieved by defining path integrals perturbatively 
as a limit D — > 1 of a D-dimensional functional integral. 
The perturbative calculations in Q were performed in 
momentum space, where Feynman integrals in a contin- 
uous number of dimensions D are known from the pre- 
scriptions of 't Hooft and M. Veltman 0]. For many 
applications, however, momentum space calculations are 
rather clumsy, for instance if one wants to find the effec- 
tive action of a field system in curvilinear coordinates, 
where the kinetic term depends on the dynamic variable. 
Then one needs rules for performing temporal integrals 
over Wick contractions of local fields. 
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II. PROBLEM WITH COORDINATE 
TRANSFORMATIONS 

Recall the origin of the difficulties with coordinate 
transformations in path integrals. Let q(r) be the eu- 
clidean coordinates of a quantum-mechanical point par- 
ticle of unit mass in a harmonic potential m 2 q 2 /2 as 
a function of the imaginary time r = —it. Under 
a coordinate transformation q(r) — > 4>(t) defined by 
q( T ) = f(<P{ T )) = <P( T )+J2^=i a n4> n (j), the kinetic term 
q 2 (r)/2 goes over into 2 (t)/' 2 (0(t))/2. If the path inte- 
gral over 4>{t) is performed perturbatively, the expansion 
terms contains temporal integrals over Wick contractions 
which, after suitable partial integrations, are products of 
the following basic correlation functions 



A(r - /) = <0(t)0(t')> = 
d T A(r - t') ee (4>(t)4>{t')) = 
d T d T ,A(r-T') ee (4>{t)4>{t')) = 



(1) 

(2) 
(3) 



The right-hand sides define the line symbols to be used 
in Feynman diagrams for the interaction terms. 
Explictly, the first correlation function (pi) reads 



A(r,r') = — e -m\r-r'\ 
2m 



The second correlation function 
1 



5 r A(r,r') 



-A-r 



(4) 

has a discontinuity 
(5) 



T '\ e -m\r-T'\ 



where 



e(r - r') = 2 f dr"5(r" - r') 

J — OO 



(6) 
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is a distribution which has a jump at r = t'. The third 
correlation function (^) contains a ^-function: 

d T d T , A(r, t') = S(r - r') - ^"H^t'I , (7) 

The temporal integrals in r-space over products of such 
distributions are undefined. 

In our previous paper we have shown that it is possi- 
ble to define unique perturbation expansions which lead 
to a reparametrization invariant theory. This is done 
by extending the path integral to a L>-dimensional func- 
tional integral, and by performing the perturbation ex- 
pansion in D-dimensional momentum space, with a limit 
D — »• 1 taken at the end. 
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In this note we want to set up an extension of distri- 
bution theory which allows us to define the same theory 
by doing the integrals over products of the correlation 
functions in r-space. This requires regularizing the cor- 
relation functions in such a way that all r-integrals yield 
the same results as the momentum integrals in the limit 
D — ► 1. The result is a reparametrization invariant per- 
turbative definition of path integrals via Feynman inte- 
grals in r-space. 



III. MODEL SYSTEM 

As in our momentum space treatment in Ref. we 
shall prove the reparametrization invariance for a simple 
but typical system, consisting of a point particle in a 
box of size d. We shall do this by performing explicit 
calculation up to three loops. 

The ground-state energy of this system is exactly 
known, E — it 2 /2d 2 . In a recent paper we have shown 
that this energy can be obtained perturbatively to any 
desired accuracy by performing perturbation expansions 
in powers of g = ir 2 /d 2 , and taking the limit g — > oo of 
this series , following the strong-coupling theory devel- 
oped in Ref. . 

In the model system, the three-loop expansion for the 
ground state energy reads 



E 



2 4 16 m 



(8) 



where m regularizes the Feynman integral in the in- 
frared. The ground state energy is obtained from the 
limit Tn — ► oo. 

After the coordinate transformation of the free-particle 
action in the path integral by q(r) — /(4>(t)) — <j) — 
gcf) 3 /3 + g 2 4> 5 /5 — ■ ■ • , we find for the energy the following 
graphical expansion 0] : 



E 



m 
~2 
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(9) 



The path integral is extended to an associated functional 
integral in a D-dimensional coordinate space x, with co- 
ordinates x^ = (t, x 2 , X3, . . .), by replacing 2 (r) in the 



d/dx^. 



kinetic term by {d^(j){x)) , where <9 P 
extended coordinate space, the correlation function 
becomes 



In the 



A(x) 



d D k 



, ikx 



(2ir) D k 2 + m 2 



(10) 



In writing down the expansion (^) we have ignored 
all diagrams coming from the Jacobian of the coordi- 
nate transformation, since these carry a prefactor 5^ (0) 



which, according to a basic rule of t'Hooft and Veltman 
, vanishes in dimensional regularization. 

In our calculations, we shall encounter generalized 5- 
functions, which are multiple derivatives of the ordinary 
(5-function: 

Si%Jx)=d^J^(x) 

= f d D k{ik)^...{ik)^e lkx , (11) 



with d Ml ... Mn =d m ... d^, and with d D k = d D k/{2ir) D . 
In dimensional regularization, these vanish at the origin 
as well: 

(0) = / d D k(ik) m . . . (ik)„ n = 0, (12) 



thus establish the contact with our previous momentum 
space discussion B. 



IV. FEYNMAN DIAGRAMS AND 
DISTRIBUTIONS 

We shall set up rules for calculating the singular ex- 
pansion terms in (^|) which will yield the following values: 



OO = -A(0)A w (0) 



D-+1 4 



V=-A^(0)A w (0) d = i -^, 



(13) 
(14) 



OOO = A 2 (0) J d D xAl v {x)=A\Q) J d D xA^(x) 



(1 + D/2)m 2 A\Q) 



d— l 16m' 



(15) 



OOO = A 2 1M (0) J d D x A 2 (x) = (1 — D/2)m 2 A 3 (0) 

(16) 



1 



d->i 16m' 



COD 030= -A(0)A^(0)Jd D xAl(x) 



= -(D/2)m 2 A 3 (0) = - 



1 



D-ti 16m 



(17) 



d D xA 2 {x)Al v {x) 



7L>)A 3 (0) 



3(3L> - 4) 
+ (Z? + 4)m 2 / d D xA i {x) 



£>-i 32m 



(18) 



d D xA{x)A^(x)A v (x)A^{x) 

2 

'(5L>-8)A 3 (0) 



6(3L> - 4) 



2 



-2(D-4)m 2 J d D xA 4 (x) 



1 



d^i 32m 



, (19) 



d D x*i(x)AUx) 



2(D-2)A 3 (0) 



3(3£> - 4) 



1 



d^i 32m 



(20) 



Summing up all contributions, we recover (g), thus 
confirming the invariance of the perturbatively defined 
path integral under coordinate transformations. 



V. BESSEL FUNCTION REPRESENTATION 

In D = 1 — e dimensions, the distributions in the 
Feynman integrals can be expressed in terms of modified 
Bessel functions K v (z). The basic correlation function in 
D-dimension is obtained from the integral in Eq. (|lo| ) as 



A{x) = c D z x - D ' 2 K^ D/2 {z), 



(21) 



where z = m |x| is reduced length of x^, with the usual 
euclidean norm \x\ — \/ x\ + . . . + x 2 D , and Ki_ D / 2 (z) 
is the modified Bessel function. The constant factor in 
front is cjj — m D ~ 2 / (2tt) d / 2 . In one dimension, the cor- 
relation function ( pl| ) reduced to (|j). The short-distance 
properties of the correlation functions is governed by the 
small- z behavior of Bessel function at origin M: 



K v (z) « -T(u)(z/2Y 

z«0 z 



Rc v > 0. 



(22) 



In the application to path integrals, we set the dimen- 
sion equal to D = 1 — e with a small positive e, whose 
limit e — > will yield the desired results in D = 1 dimen- 
sion. In this regime, Eq. (^2|) shows that the correlation 
function ( [2l|) is regular at the origin: 



A(0) 



,D-2 



(4 7 r)^/ 2r ( 1 2)0=1 2m' (23) 



The first derivative of the correlation function (21), 
which is the D-dimensional extension of time derivative 
(0), reads 



A M (z) = -c D z 1 - D ' 2 K D/2 {z)d^z 1 



(24) 



where d^z = mXfj,/\x\. By Eq. (22), this is regular at the 
origin for e > 0, such that the antisymmetry A M (— x) — 
-A^(x) makes A M (0) = 0. 

In contrast to these two correlation functions, the sec- 
ond derivative 



A^(x)=A(x) (d»z)(d lJ z) + -^-z D l 2 K D/2 {z)d^z 



„2-£> 



£)-2 



(25) 



is singular at short distance. The singularity comes from 
the second term in (Ea): 



d^z 



2-D 



= (2-D) 



m1 '" f 2.i) 



•r 



D 



which is a distribution that is ambiguous at origin, and 
defined up to the addition of a 6^ (x)-function. It is reg- 
ularized in the same way as the divergence in the Fourier 
representation Eq. (12). Contracting the indices fi and v 
in Eq. (Effl, we obtain 



(27) 



d 2 z 2 - D = (2 - D)m 2 - D S D S^(x), 



where Sd = 2n D l 2 /T(D /2) is the surface of a sphere in 
D dimension. Contracted Eq. (Ea) 



A MM (ar) -m 2 A(x)-c D m 2 ' S D \r (D/2) 2 D ' 2 5™ (x) 



= m 2 A(x) -S iD) (x), 



(28) 



coincides with the definition of the correlation function 
by the inhomogeneous field equation 



(-dl+m 2 )q(x) =SW{x) 
From Eqs. (|l|) and ^ we see that 

A MM (0) = m 2 A(0) 



m 



(29) 



(30) 



which determines the Feynman diagrams ( p_3[ ) and (|14j). 

A further relation between distributions is found from 
the derivative 

A M „ (x) = d v [-6^ (x) + m 2 A(x) 

+mS D [AOzOM^- 1 ^)] S^(x) = d„A xx (x). (31) 

Let us now turn to the calculation of Feynman integrals 
(|l5|)-(p0|) over products of distributions. 

VI. INTEGRALS OVER TWO DISTRIBUTIONS 

Consider now the integrals over products of two such 
distributions. If an integrand depends only on 

we may perform the integrals over the directions of the 
vectors 

/>OC 

d D xf(x) = S D drr^fir), r = \x\. (32) 



_ Thus we can calculate directly 



d u x A 2 {x) = mT u c 2 D S D I dzzK(_ D/2 (z) 
1 



™ Cd S d - (1 - D/2) r (1 - D/2) T (D/2) 



r D c 2 

2-D 
2m? 



A(0). 



(33) 



3 



and 



J dPx A*(s) = m 2 ~ D cl S D f 



dzzK 2 D/2 (z) 



,2-n 2 



D 



c D S D -T (1 + D/2) r (1 - D/2) = -A(0), (34) 

which together with ( |30| ) explains the values given for the 
Feynman integrals in Eqs. ( |l6|) and (|l7|). 
Note that due to the relation S: 



K D / 2 {z) = -z 



_ _ 7 Dji-\ 

dz 



(35) 



the integral over z in Eq. (|3J) can also be performed by 
parts, yielding 



d D xA 2 (x) 



7 2-D 2 



5 B (// 2 K, /2 ) (V-^ 2 i^ A 



m 



^lA^a;) = A(0) - m 2 / d^cc A 2 (x). (36) 



The upper limit on the right-hand side does not con- 
tribute because of the exponentially fast decrease of the 
Bcssel function at infinity ||. 

Using the explicit representation (H|), we calculate 
similarly the integral 

J d D xA 2 vu (x) = J d D xAl„ (37) 
= m 4 J d D x A 2 {x) - m A - D c 2 D T {D/2) T (1 - D/2) S D 
= m 4 / d D x A 2 (x) - 2m 2 A(0) = - (1 + D/2) m 2 A(0) . 



The first equality follows from two partial integrations 
and Eq. (|3l|). In the last equality we have used (33). 
We have omitted an integral containing modified Besscl 
functions 



(D-l) 



dzK D/2 (z)K 1 _ D/2 (z) + 

d r°° , - 

2 J 



dzz 1 K 2 D/2 (z) 



since this vanishes in one dimensions as follows: 



-r(i- £ /2)[r(e/2) + r(- e /2)] e 2 r( £ ) = 

4 e— >0 



(38) 



0. 



The result (|37]) explains the value stated for the Feynman 
integral in Eq. (|l5|). 

The results (p3"|), (|34|), ( |3^ ) can be used to derive a 
fundamental rule that the integral over the square of the 
^-function vanishes. Indeed, solving the inhomogenous 
field equation (pq) for 5^(x), and squaring it, we obtain 



d D x [8 D {x)Y = m A j d D x A 2 (a;) 
+2m 2 J d D x A 2 (x) + J d D x A 2 vv {x) = 0. (39) 
Thus we may formally calculate 

d D x S D (x)S D (x) = S D (0) = 0, (40) 



pretending that one of the two <5-functions is an admissi- 
ble test function of ordinary distribution theory. 



VII. INTEGRALS OVER FOUR DISTRIBUTIONS 

We now calculate the integrals over products of four 
distributions required for the diagrams (|l8|)-(f20|). These 
integrals are straightforward in D = 1 dimension, as long 
as they are unique. Only ambiguous cases require a cal- 
culation in D = 1 — e dimension, with the limit e — > 
taken at the end. 

A unique case is 



d D xA i {x) =c D m- D S 



D 



dz z 3 ' D Kf_ n/2 {z) 



Dal 

Similarly, we calculate 
J d D xA 2 (x)Al(x) 



l-D/2\ 
1 



4 - lQ fV/ 3 D \r(n\ 
Ci m oitttI x — — I (D) = 
1 2 4 V 2 2 K 1 d^i 32m 5 



(41) 



m 
1 



2-D A 



c%S D / dzz^ K D/2 (z)Kt D/2 (z) 
o 



,2-D 4 c 
1 C D JD 



2- D - 1 T(D/2)T 3 (1- D/2) 
A 3 (0) -m 2 / d D xA 4 {x) 



'■D/2 



1 



d^i 32m 3 ' 



(42) 



Using the expressions ( p4| ) and (p5j), we find for the inte- 
gral in D = 1 — e dimensions 



d D x A(x) A M {x)A v (x)A^ (x) 



m 



d u xA'(x)Af t (x)+I D , 



(43) 



where Id denotes the singular integral 

"d s d 



I D = (D-l)m i - D c D S L 



/ dzz 2 - D K x _ Dl ^{z)Kl l2 {z). 
Jo 



(44) 



In spite of the prefactor D — 1, this has a nontrivial limit 
for D — > 1, the zero being compensated by a pole from 



4 



the small- z part of the integral, at z = 0. In order to 
see this we use the integral representation of the Bessel 
function [fLlt 



1 



K u {z)=-k-^{z/2)- v T^-+v 

POO 

x / d£(coshf)~ 2,y cos(zsinhi). (45) 
Jo 

In one dimension where v — 1/2, this becomes simply 
K 1/2 {z) = ^i:/2ze- z . For v = D/2 and v = 1 - D/2 
written as v = (1 =p e)/2, it is approximately equal to 

K ilTE)/2 (z) = (z/2)-^)/2 r (l T |) (46) 

/•CO 

■e~ z ±e / dt(coshi) -1 ln(cosht) cos(zsinht) , 
Jo J 

where the t- integral is regular at z = [Q. After sub- 
stituting (|i^) into (^i|), we obtain the finite value 



It 



'2 



i-D A a \ c 



x — r (1 + e/2) T 3 (1 - e/2) x 2~ 5e T(2e) 
/ 1 



V 2m7r 2 / 8 



16m 



(47) 



The prefactor D — 1 = — e in ( |44| ) has been canceled by 
the pole in r(2e). 

The nontrivial nature of the integral Id was first ob- 
served in another form in the momentum space calcula- 
tions of Ref. (J] , where Id appeared in the integral 



d D xA\x) [A^(x) - A 2 A (*)] = -2I D , (48) 



Indeed, using the Bessel expressions (|24|) and (|2S]), we 
find 



d D x A 2 (x) [A%(x) - Al x (x)] = -(£> - l)m 



4-r> 



d: (^ 1 - D / 2 ^ 1 _ fl/2 J ^f>/ 2 , 



2 d 



(49) 



> 2 d 



(50) 



x c D S D / 
Jo 

and a partial integration 

dz ( z i-°/2 Kl _ D/2 y ^K 2 D/2 

= 2 dzz 2 - D K 1 _ D/2 K 3 D/2 
Jo 



establishes contact with the integral ( pp for Ip. 

Knowing J^, we also determine, after integrations by 
parts, the integral 

J d D xAl(x)Al(x) = -3m 2 J d D x A 2 (x)A 2 (x) - 2 I D . 

(51) 



It remains to calculate one more unproblematic inte- 
gral over four distributions: 



d D xA 2 {x)A\ x {x) 

2m 2 A 3 (0) + m 4 / d^xA 4 ^) 



(52) 



D=l 



32m 



Combining this with ( f47|) and (pl| ) we find the the Feyn- 
man diag ram (IS). The combination of ( |4^ ) and j5l| ) 
with (^) and (42), finally, yields the diagrams ( |l9|) and 
(P0|), respectively. 



VIII. SUMMARY 

In this note we have shown that with appropriate 
Bessel function representations, we can evaluate inte- 
grals over products of distributions in configuration space 
which reproduce the results of dimensional regulariza- 
tion. This ensures the invariance of perturbatively de- 
fined path integral under coordinate transformations ob- 
served in H. 
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p. 85, or Ref. M, p. 83, Formula 27. 



See Ref. [llj, Formulas 3.511.1 and 3.521.2. 
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